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ABSTRACT 

We propose a new mechanism to generate primordial curvature perturbations, based on 
the resonant production of particles during inflation. It is known that this phenomenon 
can trap the inflaton for a fraction of e-fold. This effect is governed by the mass of the 
produced particles and by their coupling to the inflaton, parameters which can depend 
on the expectation value of other fields. If one of such additional fields - a modulaton - 
is light, then its fluctuations, acquired during the earlier stages of inflation, will induce 
a spatial modulation of the trapping, and thus of the end of inflation, corresponding to a 
curvature perturbation. We calculate the power spectrum, bispectrum and trispectrum of the 
curvature perturbations generated by this mechanism, taking into account the perturbations 
due to the inflaton fluctuations as well. We find that modulated trapping could provide 
the main contribution to the observed power spectrum and lead to detectable primordial 
non-gaussianities. 
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1 Introduction 



Soon after being proposed as a mechanism to overcome the shortcomings of the standard hot 
Big Bang scenario, inflation turned out to come with a very rich bonus: it provides the seed 
inhomogeneities that eventually will evolve to form the structures we inhabit (see e.g. [1]). 
In this standard picture, there is one effective degree of freedom, the inflaton, that fulfills 
the double task of driving inflation (via its zero mode) and of generating the primordial 
spectrum of perturbations (via its quantum fluctuations). 

More recently, it has been realized that these two purposes can be fulfilled by different 
fields and one can envisage scenarios where primordial perturbations are initially stored in the 
fluctuations of a second, initially subdominant, scalar field and subsequently transferred to 
the dominant matter component of the Universe. In more technical terms, this corresponds 
to a transfer of an initial isocurvature, or entropy, mode into a final curvature, or adiabatic 
mode. This transfer can occur (i) during inflation in multi- inflaton scenarios; (ii) just at 
the end of inflation, like in the modulated reheating scenario [2], where the subdominant 
scalar field controls the decay of the inflaton into ordinary matter; (Hi) long after the end 
of inflation, like in the curvaton scenario [3], where the initially subdominant curvaton field 
comes to dominate after reheating as it redshifts more slowly than radiation. 

Another aspect of inflation which has been actively investigated in the last few years is 
the resonant production of particles, which can arise due to the coupling of the inflaton to 
other fields. Indeed, if the inflaton is coupled to a field x? either bosonic or fermionic, then 
X particles can be produced by resonant effects whenever their effective mass crosses zero as 
the inflaton evolves. The most studied example is preheating [Ij (see also [5]), where particle 
production occurs while the inflaton is oscillating at the bottom of its potential. However, 
this effect can also take place during inflation as pointed out in [6]. In the latter case, 
the backreaction of the produced particles induces a slow down of the inflaton 0, effectively 
trapping it for a fraction of e-fold before the quanta of x are diluted away by the expansion of 
the Universe. One consequence of this slow-down of the inflaton is the generation of features 
in the inflationary spectrum O [71 [8]. The same mechanism can also help halt moduli at 
points with enhanced symmetry [9]. More recently, it has been shown in [10] that repeated 
trapping events can slow down the inflaton enough to lead to slow-roll inflation even in the 
presence of a steep inflationary potential. 

In the present paper, we consider a scenario where the intensity of the trapping depends 
on an additional light scalar field, a, which we will call a 'modulaton' (this field does not 
contribute significantly to the energy density, neither during inflation, like in multi-field 
inflation, or after inflation, like in the curvaton case). More precisely, we will assume that the 
mass of X, or its coupling to the inflaton, depends on the modulaton cr, which is light during 
inflation and thus acquires a quasi scale-invariant spectrum of super-Hubble fluctuations. 
Consequently, in super-Hubble regions of the Universe with different values of a, inflaton 
trapping will occur at different times and/or will be more or less strong. The duration of 
inflation will thus vary from one region to another. This way the perturbation in the x 
modes will be transferred into a curvature perturbation. 

This idea of using a modulaton is not new. This is indeed the key ingredient to the inho- 
mogeneous or modulated reheating scenario [2j , where the decay rate of the inflaton depends 
on the modulaton cr. It was extended recently [11] to modulated preheating, where a phase 
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of preheating takes place instead of perturbative reheating. In this case, the modulaton en- 
ters into the couphng g between the inflaton and the produced particles and the fluctuations 
of the modulaton are transferred into curvature perturbations, because both the duration 
of the preheating phase and the density energy of the scalar field at the end of preheating 
depend on g (similar ideas are discussed in [T2]). 

While using the ideas of modulaton and particle production, our scenario is very different 
from the previous ones because the crucial effect here is the backreaction of the particle 
production on the motion of the inflaton or, more generally, of any scalar field dominating 
the energy density in the Universe. After they have temporarily trapped the inflaton, the 
produced particles become cosmologically irrelevant as they are rapidly diluted by expansion. 

Using the 6N formalism, we compute the amplitude of the curvature perturbations gener- 
ated by modulated trapping and we compare their contribution in the power spectrum with 
the contribution due to the usual inflaton perturbations. We distinguish the case where the 
coupling depends directly on the modulaton from the case where the modulaton affects only 
the critical value of the inflaton at which particle production takes place. In the latter case, 
the contribution from modulated trapping is usually small in the simplest models of inflation 
but it can become dominant if the Hubble parameter at particle production is much smaller 
than when the modulaton fluctuations were generated. In the former case, modulated trap- 
ping can easily dominate. We also investigate the primordial non-Gaussianities generated 
by this mechanism. More precisely, we compute the bispectrum and the trispectrum and we 
find that modulated trapping can lead to a level of nongaussianities that would be detectable 
with future experiments. 

The plan of our paper is as follows. In section 2 we present the model and review the 
main results concerning the impact of resonant particle production on the evolution of the 
slow-rolling inflaton. In section 3 we derive the linear perturbations generated by modulated 
trapping and compare their contribution with that due to the inflaton fluctuations, for 
different cases. In section 4 we investigate non-Gaussianities by computing the bispectrum 
and trispectrum. Finally, we conclude in section 5. 

2 Model and homogeneous solution 

We consider a model where the inflaton is coupled to other fields Xy which can be either 
bosonic or fermionic. In addition to the usual kinetic terms and the self-interaction potential 
V {(f)) of the inflaton, we assume that the Lagrangian contains, in the bosonic case, a term 
of the form 



In both cases. A/" denotes the number of species of particles with the same mass m and the 
same coupling A. Taking into account the coupling to the inflaton, these particles have an 
effective mass 





and, in the fermionic case. 




(2) 



M(0) 



m — A0, 



(3) 
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so that there is a critical value of 0, denoted by 

0^ = m/A (4) 

where M (0) vanishes. As discussed in |1] , resonant particle production will occur when cf) 
crosses 0*. In this section we will review the main results of [6j, where resonant production (of 
fermions) during inflation was first studied. As in [6j, we will focus on fermion production, 
keeping in mind that the scenario where scalars are resonantly produced gives the same 
results. 

When M (0) crosses zero, the particle occupation number suddenly increases from zero 
to the value 

_X3/2 

where we use the subscript * to denote the time of particle production. The number of 
particles is then diluted by the expansion so that 

n{t) = (—\ ^e{t-Q , (6) 

where is the Heaviside distribution. The backreaction of the particle production can be 
estimated by using the Hartree approximation in the equation of motion of the infiaton, 
which after substituting the particle number obtained above, becomes 

4> + 3H<p + V'{<j)) =ArX{xx) = >^nJ-) ^Q{t-Q (7) 

Since the production and subsequent dilution of the particles occurs during a fraction of 
e-fold, we will assume H = if*, V {(p) = l^' (0*) to be constant during the entire process. 
One can then easily integrate the equation of motion for 0. Denoting 

A0(t) = 0(t, A ^ 0) -0(t, A = 0) (8) 

the difference between the solutions with and without particle production, one finds 

A0(t >t,) =ArAn,e-=^^*(*-**) (t-t,) . (9) 

We have assumed here implicitly that < and the positive A0 thus corresponds to a 
decrease of the absolute value of the infiaton velocity. Since particle production occurs at 
the expense of the kinetic energy of the infiaton, the sign of A0 will always be opposite to 
that of 0, so that |0 + A0| < |0|. 

This slow-down of the infiaton, illustrated in Fig. 1, leads to the generation of features 
in the inflationary spectrum since its amplitude depends on 0. This has first been pointed 
out in [6j and subsequently studied in [7] and [8]. A scenario where the resonant production 
of particles during inflation leads to features in the observable curvature spectrum has also 
been discussed in [13] . In the present work, however, we will assume that particle production 
occurs much later than when the scales of cosmological interest crossed out the Hubble radius. 
Consequently, the features in the spectrum will affect scales which are much smaller than 
those corresponding to the present cosmological window. This also applies to the fluctuations 
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Figure 1: Evolution of \(f)\ as a function of the number of e- folds for a model of chaotic 
inflation with trapping (lower curve) and without trapping (upper curve). 



generated by the re-scattering of the produced particles with the homogeneous scalar field, 
which have been recently studied in [T3]. 

Of course, the energy injected into the produced particles cannot exceed the kinetic 
energy of the inflaton before production, and the value of Acf) is limited by the fact that 
the rolling of (f) will at most be halted by trapping. In other words Max|A</)| cannot exceed 

= The behavior of the system when backreaction effects are important has been 

studied numerically in [B], where it was found that the evolution of A(/) has still the functional 
form of ([9]), but with a different overall normalization that ensures energy conservation. We 
can account for a such a behavior by defining the parameter 



(3 = Max(A0)/t;* 



^^\ 



(10) 



which must satisfy /5 1. The above expression can be used to trade A for /?. It is also useful 
to note that, if the inflaton is in slow- roll, its velocity f* just before particle production is 
given by 

v^ = ^/2e^HMp, (11) 
where e is the first of the usual slow roll parameters 




V 

and Mp = (SttG)"^/^ is the reduced Planck mass. 

The inflaton is only momentarily slowed down, as the term on the right hand side of 
quickly goes to zero after the production has occurred {t — t^, 3> 1/H^). However, the inflaton 
will have accumulated a delay with respect to its evolution without particle production: 



A0 



A0cit 



MX 



3/2 



dm 



(13) 



Consequently, inflation with particle production will end later than inflation without particle 
production, as one can see clearly on Figure [2l 
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Figure 2: Evolution of as a function of the number of e- folds with trapping (upper 
curve) and without trapping (lower curve), using the same model and initial conditions as 
in Fig. 1. 



3 Fluctuations 

We now assume that the fermions (or bosons) are also coupled to another light scalar field 
cr, which we will name the modulaton field to distinguish it from the infiaton, so that the 
effective mass 

moff = m(cr) - A(a) </> (14) 

depends on a. One possibility, which we will consider later as a specific example, is that the 
coupling to cr arises from a Yukawa coupling —g (JXXi which implies m(cr) = ga. 

The critical value for where particle production occurs now depends on a and is given 

by 

, , mia) 
X[a) 

Since the modulaton a is assumed to be light, i.e. its mass is much smaller than the Hubble 
parameter during inflation, it acquires some fluctuations on super-Hubble scales, via am- 
plification of its quantum fluctuations. As a consequence, the value of the scalar field a 
fluctuates from one patch of the Universe to another, and therefore particle production is 
slightly different in each of these patches. The time delay in the evolution of the infiaton, 
calculated in the previous section, will thus fluctuate from one patch to the other, which can 
be reinterpreted as fluctuations in the final curvature perturbation at the end of inflation. 

In order to quantify the curvature perturbations generated by the fluctuations of a, it is 
convenient to use the so-called 5N formalism |15j, based on the local number of e- folds, or 
integrated expansion, between some initial and final hypersurfaces 

N{x) = / H{t,x)dt. (16) 



In particular, the curvature perturbation on a uniform energy density final hypersurface, 
which we will denote (, is directly related to the perturbation of the number of e-folds 
defined with respect to some initial flat hypersurface, 

C = 6N = N{x) - N, (17) 
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where N is the number of e-folds in the homogeneous background spacetime. This is true 
not only at hnear order but also at non-linear order [T71 [T5] . which will enable us to 
compute the non-Gaussianities in the next section. 
During slow-roll inflation, the number of e-folds is 

'f'f H , I f^f V 



where we have used the slow-roll equations = —V'I{?>H) and 3if^ = V/Mp to obtain the 
second equality. In our case, the above expression is not valid during the bursts of fermion 
production, where the scalar field deviates from the slow-roll solution. This induces a net 
change in the number of e-folds, 

A/V = -i/.^ = ^ = ^!^ (19) 

where we have used f|T3l) . This shift in the number of e-folds depends on the scalar field a: 
either directly in A, or indirectly via the time of particle production, which depends on 
the critical value 0*(cr). 



3.1 The power spectrum 

Since the masses of both the infiaton and (by assumption) the modulaton a are smaller 
than the Hubble parameter during inflation, both and a will carry perturbations on super- 
Hubble scales, characterized by the power spectra 

Vs^ = Vs.= {^)\ (20) 



TT 



where the index k denotes the Hubble crossing {k = aHk), for the cosmological scales of 
interest . 

By expanding the number of e-folds at first order as a function of the perturbations of 
the infiaton, 50, and of the modulaton, Sa, one can write the curvature perturbation ( as 

( = SN = (5iV3io„_,oii + SAN = 50 + AiV,. Sa . (21) 

Mp Vf. 

This implies, since the perturbations in and in a are uncorrelated, that the curvature 
power spectrum is given by the sum of two contributions. 



Mp Vf, 



+ {an^Ym^p 



' (22) 



27r Mp 



where we have used (1121) and (1201) in the last equality. 

The importance of the trapping effect with respect to the usual infiaton fluctuation term 
can be quantified by its relative contribution in the total power spectrum, defined by 

^ _ ^trap ,r)o\ 

p ' (^^) 
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so that _ 

2efc(AiV,.)2M|.. (24) 



^trap 



1 - S 

One recovers the standard inflation resuh when S ^ 1. The opposite hmit is S = 1 where the 
trapping dominates. The intermediate values of S correspond to situations where both terms 
can significantly contribute. This is quite similar to mixed inflaton and curvaton models, 
investigated in [19], where both the inflaton and the curvaton fluctuations contribute to the 
observed spectrum. 

The spectral index is also affected by the fact that the two contributions have different 
scale dependence, since the inflaton contribution is proportional to -ff|/efc while the trapping 
contribution is simply proportional to One finds that the spectral index is given, in 
terms of the slow-roll parameters and of the trapping fraction H, by the expression 

- 1 = ThTfc " -^)Vk + (4S - 6) efc, (25) 

at leading order in the slow-roll parameters. It can be checked that this result is exactly the 
same as in the mixed inflaton-curvation scenario [19]. If the inflaton contribution dominates, 
one recovers the standard result — 1 = 2rik — 6ek while one finds — 1 = —2ek if modulated 
trapping dominates, like in the pure curvaton scenario. 

Let us now compute (AA^^^-)^ explicitly. Substituting the expression of v^: in terms of e* 
and into (|T9l) yields 

where e* = e(0*) and = H{(j)^) depend on a via 0*(cr) = m/X. Using 




V2 . .^ H^^ V V2t 

(r/ - 2e) , — = — = — — (27) 



e V ^' ^ ) Mpy/I " " H 2V 2M, 

and 



da da \ X 



(28) 



^"^■^ lT + 7lf^l/-(^*-3e.)(^) , (29) 



one finds 

'an 2 X ' AMp 

where the primes in the above formula denote a derivative with respect to a. We thus finally 
obtain for the power spectrum due to modulated trapping 




o AT^A -5Mp^/2^,XX' + {r],-3e,) (mA'-Am') 



271 J 1296 7:^ MpH^^/2^ V27r/ 

(30) 

In the next two subsections, we will consider two particular cases of the above general 
formula. In both cases, we will assume that the mass of the fermions depends linearly on a, 
i.e. m = ga where g is constant, as would result from a standard Yukawa coupling between 
a and the fermions. We will moreover assume that the coupling A is independent of a in the 
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first case, whereas it is of tlie form A = o jM^ such as would arise from a non-renormahzable 
term, in the second case. In the latter case, A and m have the same hnear dependence on a 
and therefore 0,, becomes independent of a . Of course, more general cases, combining both 
direct and indirect dependence on a can be envisaged, and the purpose of our two examples 
is simply to analyze separately the two types of dependence on a. 

3.2 Modulaton- independent coupling 

We now assume 

A = Ao, m = ga, (31) 

where Aq and g are constants. The trapping effect then depends indirectly on a via the 
critical value 0*(cr). 

Then, the expression fl30|) reduces to 

Kap - (^^J 9 M Ao ^ H^Mp- 

Reexpressing Aq in terms of the braking parameter jS, via ffTOj) . the corresponding amplitude 
can be written in the form 

rli = A.f:,"'{0'\ (33) 
where we have introduced the fraction /^.^ = H^/Hk and 

/ p3 \ 1/5 

A*= (y^T^J ^7Ar2/5/5'/'s. -4.5x10-3(7 Ar2/5/33/5s^, s.{e,,r],) = t-,^"' \r], - ?,e^ . 

(34) 

Let us now see whether the trapping contribution can dominate the power spectrum, i.e. 
S ^ 1. If this is the case, then the trapping amplitude (|32l) must agree with the observed 
amplitude of the primordial spectrum of perturbations 

^ 4.9 X 10"^ . (35) 

This implies the condition 

^^4.1xlO-«A;VV.f . (36) 

Moreover, we must check that the infiaton contribution is much smaller than the trapping 
contribution, which implies the condition 



- = 87^\Al » 1 . (37) 



1-S ^" ^"'^^ \Mp) 
Replacing A* by ( IMl) and Hk ( l36l) thus yields the condition 

^5/2_^^3/2^-V2^^^5/2>>g4 (38) 

on the various parameters of the model, for a domination of the trapping effect. 
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If one takes, for instance, g = 1, A/" = 100 (the same choice of parameters was considered 
in [6J) with Aq large enough to make backreaction important {(3 ~ 0.8), one finds that 
g^l'^ J\f jj'^l'^ ~ 71. Since the slow-roll parameter must be much smaller than one, this 
means that the condition fl38l) cannot be satisfied unless <C t]^ or if the Hubble parameter 
varies very much between the time when the perturbations are generated and the time when 
trapping occurs, so that f^k ^ 1- The former case, i.e. ^ r/*, can arise if inflation at 
the time of particle production behaves like in hybrid inflation. The latter case can occur 
if inflation proceeds through different phases with a hierarchy between the various energy 
scales (like in double inflation), or even if the Universe underwent several distinct phases of 
inflation, as discussed for instance in [20] or [21] . 

Even if the trapping is not the main contribution, it can nevertheless represent a non- 
negligible fraction of the overall fluctuation spectrum, when the parameters and s-t are 
not too small. For instance, in the case where e* ~ r^* ^ 0.5, ~ 0.01 and ~ 0.1, one 
obtains S ~ 0.1. 



3.3 Modulaton-dependent coupling 

We now assume ^ 

A = ^, rn = ga, (39) 

where g is constant and where M is some mass scale. In this case, AA^ depends on a 
directly via the coupling A, while the critical value 0* of the inflaton at which trapping 
occurs is independent of a. Substituting (139|) in (l30l) . we obtain 

p.„„.(_A^)V^A»^|i|?. (40) 

By comparing this expression with the analogous one of the previous subsection, i.e. (l32l) . 
we see that the effect of modulated trapping can be much stronger in the present case. The 
numerical factor 25, the different dependence on the slow roll parameters, and especially the 
fact that Mp/M can be much larger than one, all concur to increase the contribution of the 
trapping and therefore push S towards unity. 

As in the previous subsection, we can write the amplitude of perturbations due to trap- 
ping in the form 

p^i = A./;."=(^)'". (41) 

with now 

A* = 5v^(y^)'^' P'^' el/'' ^ ^ 3 X 10-' ^^'/' el^^' ^ . (42) 

This gives for example ~ 0.2 /^/^^ M^/^ hI'^/M when we take = 0.5, M = 100 and 
P = 0.8, as in the previous subsection. 

The contribution due to the trapping dominates the observed spectrum if the condition 
(137|) is satisfied. Using the present expression for A*, this condition can be rewritten as 

A/'/5^/V;//^.ei/^(^)''''»0.5, (43) 
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which is rather easily reahzed, either with a large number of species or a small ratio M/Mp. 



4 Non-Gaussianities 

Let us now study the non-Gaussianities in the modulated trapping scenario. In standard 
single field infiation, primordial non-Gaussianities are negligible but in many other scenarios, 
they could be significant. An easy way to compute the non-Gaussianities of the curvature 
perturbation is to use the Taylor expansion of the local number of e-folds in terms of the 
scalar field fiuctuations [22]. For a multi-field system, this gives 

^ = m = Ni 5/ + ]-Nij 5/ 5^^ + \nijk 5^-^ 5<^^ + . . . , (44) 

where we use the implicit summation convention for the field indices I, J,K, . . . and the 
notation Nj = dN/d(f' , Njj = d'^N/dif'dif' , etc. In Fourier space, this implies that the 
three-point function is given by 

(Ck.Ck.Cka) = i2n)''S^'\Y.^,)B^ik,MM) = NrNjNKiSifiiMM,) + 

i 

+]^NiNjNKL{5vl,5^i,{5^^ ^ 5^%,) + perms, (45) 

where the symbol * denotes a convolution product. 

If the scalar field fiuctuations are independent and Gaussian, which is expected to be a 
good approximation for infiation with standard kinetic terms and in the slow- roll limit, one 
can ignore the three-point correlation functions of the scalar fields and take into account 
only the two-point functions 

= {27^fSuS^'\\^i + k,)P(fcO, P{k) ^ '^nk), V{k) ^ 0, (46) 



The bispectrum (1451) can then be written as 

i?c(ki, k2, kg) = ^/nl [Pik^)P{k2) + P{k2)Piks) + Piks)P{ki)] . (47) 



with the non-linearity parameter 



6 _NjNjN^ ... 
5^^^ - {NKN^r 



The present observational constraints on the non-linearity parameter /nl, based on the 
WMAP 5yr data, are |23i i24j 



- 4 < f°l < 80 (95%CL) (49) 

for the local type of non-Gaussianity considered here. The Planck satellite, which has just 
been launched, is expected to reach /nl ~ 5. 
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Using the Taylor expansion up to third order, one can compute in a similar way the 
trispectrum, i.e. the Fourier transform of the connected four-point function defined by 

(Ck.Ck.CkaCkJc = {2nfs^'\Y. k,)r^(ki, k2, k3, k4) . (50) 

i 

With the same assumptions as above, the trispectrum can be written in the form [25] 

54 

T^{k^, k2, kg, k4) = tnl [P{h3)P{ks)P{h) + 11 perms]+— [P{k2)P{h)P{h) + 3 perms] 

(51) 

with 

""^^ " (NlN^Y ' " 54 (NlN^Y ^^^^ 

and where fci3 = |ki + ksj. 

After this general introduction, let us consider our particular model where the number of 
e-folds contains two separate contributions, so that the Taylor expansion up to third order 
is given by 

C = 6N= '^^^"' s<P + ... + ANJa + l-AN^Ja' + ]:AN.^Ja' (53) 

where we have ignored the second and third derivatives with respect to the inflaton, which 
give negligible non-Gaussianities. According to (j48l) . the corresponding non-linearity param- 
eter is given by 

5^""" {Nf + {AN^)f [ n ) i^N^y ^ ^ ^ 

Similarly, the coefficients of the trispectrum flSTl) are 

_(AiV)!^3_ 36 2 „ _25AiV^^3 

"^^^ " {AN^y " " 25H^^L ' " 54 {AN^^ " ^ ^ 



So far, our expressions fl54ti55p are quite similar to those obtained in the context of mixed 
inflaton and curvaton models |26] or modulated reheating [27] . 
Using now our explicit expressions ( !26l) and ( !29l) . one flnds 



AN.. 



(AiV,,)2 e/5 (5MpV2^A'+ (r/, -3e,) A0t 

X '(2Ul-8r],e,-r]l + 2V2r,Mpr]')X^(j)',^ + We,(3X'^ + 2XX") Mj,+ 
+ 2 V2r, (r/, - 3e,) (A^ 0'; + 5 A A' 0',) Mp] , (56) 

which we have expressed in terms of 0=„(o") (rather than in terms of m(cr)) for compactness, 
and where we have used eq. ffTOl) . Below, we specialize this expression for the two particular 
cases which we have discussed in the previous section. 
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4.1 Modulaton- independent coupling 



In the case A = Aq, m = go, we find from ( l54l) and ( l56l) 

, _ 5 ^ _ 2l6^-8r7,6, + 2Mpr/:v/2^-r/g 

(r^.-3 6.)2 • ^^^^ 

Tlie term B^, in tlie above formula is typically of order one, since both the numerator and the 
denominator of the above formula are proportional to the square of the slow-roll parameters. 
For instance, chaotic inflation corresponds to 77 = e, MpT]' = —r]\/2^, and thus = 2. 
Even if trapping dominates the spectrum, i.e. S ~ 1, the non-linearity parameter cannot be 
much bigger than unity. Note also that, if 77* 3> e*, one finds i3* = — 1 and thus a small but 
negative /nl- 



4.2 Modulaton-dependent coupling 

In the case A = a/M and m = ga, (i.e. 0'^ = 0, A" = 0), we find the remarkably simple 
result 

fc^^H^. (58) 

The parameter fj^^ in this case depends only on the braking parameter f3 and on the relative 
contribution of the trapping in the power spectrum. If the trapping dominates, a small 
braking index P can lead to a significant /nl: for example, P = 10~^ gives /nl — 55. 
However, one must be aware that a small (3 tends to reduce S as well. In the small (3 limit, 
S is proportional to and /nl thus scales like /?^/^. 

One can easily generalize the above result for an arbitrary coupling and one finds 



,AA"^ 



/nl = — (3 + 2^ (59) 



One can also go beyond the bispectrum and study the trispectrum. Substituting the 
third derivative of AA^ with respect to a in the second expression of (l55l) yields 



5'NL 



2e2/52 



■ AA" A^A'" 



(60) 



It is interesting to observe that, if the coupling is of the form A(cr) = (cr/M)^, then one gets 

^^^ = 2|5('-^)' "^^ = 2ll^('-^) ' = el:^^ - ^) - ^) • ^''^ 

In this case there is a simple relation between the two coefficients tnl and g^i^, which depends 
only on the power p, 

54 (5p - 2) 
^^^=25(5p-4)^^^' 

while S can be determined from tnl and /nl since S = (36/25)/NL/rNL- 
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5 Conclusions 



We have shown that modulated trapping, due to resonant particle production, provides a 
new mechanism for the conversion of isocurvature into adiabatic perturbations. Although 
we have discussed mainly the trapping of the inflaton, it is worth emphasizing that the 
same mechanism would apply to any rolling scalar field whose energy density contributes 
significantly to the total matter budget in the Universe. 

In the present work, we have computed the perturbations generated by the trapping 
effect and compared their contribution in the power spectrum with respect to the usual 
contribution from the inflaton fluctuations. We have also computed the non-Gaussianities 
generated by the trapping effect. We have found that the non-linearity parameter /nl of 
the bispectrum is proportional to the square of the fraction S of the power spectrum due 
to trapping, while the parameters tnl and (^nl of the trispectrum are proportional to S^. 
This means that the non-Gaussianities are suppressed if S is small. To explore whether 
modulated trapping could account for most of the observed perturbations, it is convenient 
to distinguish two main cases. 

In the first case, where the modulaton field affects only indirectly the particle production, 
the trapping contribution in the observed spectrum of fluctuations is expected to be small, 
except in two types of scenarios: if trapping occurs at a Hubble scale which is much smaller 
than the Hubble scale when the perturbations were produced, i.e. H^: <C Hk, or if trapping 
occurs in a slow-roll phase where the slope of the potential is much smaller than its curvature, 
i.e. e* <ti 1]^. Moreover, in this first case, the parameter /nl is at most of order unity (if 



More interesting is the second case where the coupling A depends directly on the mod- 
ulaton. It is then easy to obtain a domination of the trapping contribution in the observed 
signal. Moreover, the non-linearity parameters /nl, tnl and (y^NL, are given by very simple 
expressions that are inversely proportional to the braking parameter (3 or to its square. 

To conclude, modulated trapping can be an efficient mechanism to generate primordial 
cosmological perturbations, with the possibility to produce significant non-Gaussianities that 
could be detected in the near future in the CMB observations. It would be worthwhile to 
study further this mechanism in specific models of inflation embedded in realistic particle 
physics setups. 
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